Spin dynamics in the stripe phase of the cuprates 
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Within a model that supports stripe spin and charge order coexisting with a dj;2_j,2-wave super- 
conducting phase, we study the self-consistently obtained electronic structure and the associated 
transverse dynamical spin susceptibility. In the coexisting phase of superconducting and static stripe 
order, the resulting particle-hole continuum can strongly damp parts of the low-energy spin wave 
branches. This provides insight into recent inelastic neutron scattering data revealing the dispersion 
of the low-energy collective magnetic modes of lanthanum based cuprate superconductors. 

PACS numbers: 74.72.-h, 74.25.Ha, 74.25.Jb, 75.40.Gb 



The electronic properties of the underdoped cuprates is 
dominated by competing instabilities and coexistence of 
several ordered states. An example is given in the accu- 
mulating evidence the doped holes tend to self-organize 
into one-dimensional rivers of chargeQ. 

For instance, in Nd-doped La2-a;Sra;Cu04 neutron 
scattering (NS) revealed a quartet of incommensurate 
(IC) elastic magnetic peaks at ((1 ± 8)it, tt), (tt, (1 ± 5)-^) 
and Bragg charge peaks at (±27r5,0), (0, ±27r(5)l3|. This 
is consistent with the doped holes forming ID domain 
walls oriented along the Cu-0 bonds in the Cu02 planes. 
These stripes are separated by l/i5 (in units of the lat- 
tice spacing a) and the staggered magnetization gains an 
extra phase shift tt when crossing a stripe. Later, it was 
realized that d-wave superconductivity (dSC) coexists 
with the static stripe order 0. The NS results from pure 
La2-a;Sra;Cu04 exhibits similar IC peakspi|. For these 
materials, Bragg peaks are observed for 0.02 < a: < 0.13, 
whereas for a; > 0.13 a small doping dependent spin gap 
opens in the magnetic excitation spectrumQ. 

The dispersion of the IC peaks, i.e. (5(w), between 0-40 
meV in optimally and underdoped LSCO was measured 
recently by Christensen et ai. 6]. It was found that both 
in the pseudogap and dSC state, the IC peaks disperse to- 
ward (tt, tt) as the energy increases but with no sign of an 
intense resonance feature at (tTjTt). Furthermore, the IC 
peaks broaden as the energy increases. A similar disper- 
sion was found in La2_a;Baa;Cu04 by Tranquada et al. 7\. 
In La2_2:Baa;Cu04 and YBa2Cu306+K it was further 
found that above the (tt, tt) crossing, the spin response is 
dominated by four peaks rotated 7r/4 relative to the low- 
energy IC peaks 0,0. These high-energy peaks disperse 
to larger wave vectors with increasing energy transfer. 
In YBCO the low-energy spin response is dominated by 
the commensurate (C) (tt, tt) resonance which disperses 
downward with decreasing energy This points to a 
degree of ubiquity in the spin fluctuation spectrum of 
the cuprates. The main difference between the materials 
appears to be the size of the doping dependent spin gap 
and the intensity distribution along the spin branches. 
These experiments have sparked new theoretical studies 
dealing mostly with the high-energy response [lfij|. 



In this Letter we report the self-consistent results of 
the electronic structure and the corresponding transverse 
spin susceptibility within a model that supports static 
IC spin and charge density wave solutions. In particular, 
we focus on the low-energy spin dynamics and the influ- 
ence of the dSC on the intensity distribution of the spin 
branches when it coexists with static spin and charge or- 
der. These studies are motivated by the above-mentioned 
new experimental insight, and the strong evidence for 
stripes in the lanthanum based materials. A similar ap- 
proach was used to study the phonon anomalies caused 
by collective modes in the stripe phase of the 2D Hubbard 
model without allowing for dSC order [ll). 

Several previous studies of magnetic modes in d-wave 
superconductors have started from spatially homoge- 
neous phases m [13. In contrast, recent spin-only models 
have also been proposed to describe the spin dynamics 
in the IC stripe phase^^l- Here, we bridge these two 
apparently complementary approaches by solving self- 
consistently the following minimal model defined on a 
2D lattice 

H = ~ {^^3^\>yC]<y + H.C.) + ^ (C/ {flia) - m) flic^ 

+ E(^'^4t4+H-c.). (1) 

Here, creates an electron of spin a on site i, tij ^ t, t' 
denote the first and second nearest neighbor hopping in- 
tegrals, /i is the chemical potential, and fiic^ — cl^Cia is 
the occupation number on site i. The model ^ is the 
mean-field version of the extended Hubbard model with 
onsite repulsion U and nearest neighbor attraction V and 
is aimed to mimic essential features of phases of coexist- 
ing spin, charge and dSC order. The nearest neighbor 
attraction V triggers the singlet dSC at the mean-field 
level, Aij = V {{ci^Cji) — (ciiCj^)). This approach has 
previously been used extensively to gain insight into the 
electronic structure in phases of coexisting order 'l3,'l6|. 

As is well-known, away from half-filling the Hamilto- 
nian produces inhomogeneous spin and charge order 
= hii^i^) ^ (^a)) and Pi = {{fill) + {hii)) with the 
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modulation period of being exactly twice the period 
of Pi . This is the stripe phase at the mean-field level. In 
general, i.e. when t' / 0, the stripes are metallic since 
the mid-gap states cross the Fermi level [lij. 

For an ordered array of stripes with spin periodic- 
ity N we can divide the real-space lattice into smaller 
supercells of size N. Then, the Hamiltonian can be 
diagonalized by the Bogoliubov transformation, — 

Enki<kA^^)e~"^'^'lLa+'^'^nkwir^)e'^^'^'%l,a), where 

Tj denotes a site within the supercell which in turn is 
positioned at R^. The wave vectors k belong to the cor- 
responding reduced Brillouin zone, and a = +(— )1 for 
up(down) spin. This results in a set of Bogoliubov-de 
Gennes equations to be diagonalized for each k. The 
self-consistency is enforced through iteration of the rela- 
tions, 1 — nil = '^i^{fiia) , nh is the hole doping, and 
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^ij^Yl Kk^(*)""k<T(j)/(^nka) (3) 
71 k 

- U«k<T(«Xk<T(i)/(--E'nko-)] • 

As usual, f{E) = [1 + exp{Ef3)]^^ denotes the Fermi 
distribution function with f3 — 1/kT. 

Below, we report the spin dynamics obtained from the 
stable configurations both with and without dSC when 
U = 4.Qt, Ufi = 0.125 and N = 8. In principle, for a given 
set of parameters C/, V, t' , and rih the system may prefer 
a spin period different from N = 8. By varying rih we 
have checked that this possibility does not qualitatively 
alter the results below. 

A typical example of the self-consistent results for the 
spin (charge) density S'f (pi) and the pairing poten- 
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shown in Fig. Clearly, this shows the expected anti- 
phase stripe ordering with concomitant modulations of 
A^. The bond-centered solutions are found to have en- 
ergies slightly lower than the site-centered stripes. 

Turning briefiy to the electronic structure obtained 
from these self-consistent solutions, it is well-known that 
for the spectral weight, /(k) = ^(k, Li;)c?w, the 

stripes generate weight around the antinodal regions [iH 
Il7| . Here yl(k,w) is the single-particle spectral function 



and Aw is an integration window below the chemical po- 
tential p. For instance, in Fig. ^ we show /(k) for 
< kx, ky < IT for the same parameters used in Fig. ^ 
but allowing for stripe disorder 18J. This figure is strik- 
ingly similar to the recent ARPES data by Zhou et al. ^ 
Furthermore, the stripe ordering causes a modulation of 
the pairing potential as seen in Fig. ^ which influences 
the Fourier transform of the LDOS A^q(w)|i3,[i^]. For in- 
stance, as seen from Fig. the real part of TVq. {uj) ver- 
sus energy lu at the charge ordering vector Q* — (^,0) 
exhibits two low-energy zero crossings in agreement with 
recent STM measurements by Howald et aJ.ET 




FIG. 1: (a) Bond-centered spin density { — lY'^Sf (red), charge 
density pi (blue) and pairing potential Af (green) versus site 
obtained when U = 4.0t, V = 2.0t, t' = -O.S7t, (b) the spec- 
tral weight /(k) (Acl) — O.lt) versus k for disordered stripes, 
(c) plot of the real part of A^'q* (lo) versus energy uj. 



Motivated by the agreement of the self-consistent 
mean-field solutions of Eqn. and the mentioned elec- 
tronic probes, we turn now to our main topic: the collec- 
tive spin dynamics in the stripe phase of Eqn. The 
stripes explicitly break the SU(2) spin rotation symme- 
try and hence associated Goldstone modes are expected 
in the transverse spin susceptibility (pt?^)- Below, 
we investigate the effect of the charges and the dSC or- 
der on the dynamic part of the spin modes. An explicit 
calculation of the Fourier transform of x'' (r^ , fj , t) = 
— {TrSj'{T)S^{0)) shows that to Gaussian order the re- 
sponse is given by the diagonal elements of the N x N 
matrix 

X+'(q,c^) - xJ-(q,o;) (1 - Uxi-{q,uj)y' . (4) 

Here the matrix elements of the bare susceptibility 
Xo ~ (q> ^)qQ' are given by 



Xo (q,t^)QQ' 



-y 
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, N 1 - /(-E^riko-) - /(-Emk-l-qcr) , , , /(-^nko-) + /(^'mk-l-qtr) - 1 /^.x 



- Em'k+qa + £^nk(T + iF 



— Emk+qcr — Enkcr + iF 



fjEnka) — /(-Emk+qj) 



f{Emk+qa) — f (Enka) 



„iq-(rj— ri)+iQ-rj — iQ'-r; 



where Q are N reciprocal lattice vectors of the supercell lattice. The coefficients ai{u,v) and bi{u,v) are given by 
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the following combinations of the coherence factors u and v 

ai{u,v) = <ka(r»)"mk+qa(rO (w„k^(j>™k+q^(rj) -u„k<T(rj)<„k+q„(i"i)) , (6) 
hi{u,v) = <k^(rj)w„k+qo^(rO (unk<T(rj)<,k+q5^(ri) + w«ka(rj)<k+qa(ri)) ■ (7) 



These coefHcients are independent of the reciprocal vec- 
tors Q since the eigensystem of the Hamiltonian for a 
given k is invariant under shifts k — > k ± Q. The re- 
maining factors, a2{u,v) and b2{u,v), are obtained from 
ai(u, v) and v) by interchanging Ua and <-> 

Uct, respectively. In Eqn. ((SJ, the real-space sum Vi^Vj, 
extends over a single supercell of size N whereas k be- 
longs to the reduced BZ. Below, we set kT = T = O.OOSt. 
Due to the large number of diagonalizations involved in 
the sum in Eqn. ((SJ, we are unfortunately unable to ob- 
tain high-resolved 2D constant energy cuts through the 
BZ. However, as shown below, the explicit to dependence 
in Eqn. lO allows us to calculate e.g. {qx,uj) plots. 

Without spin order {U = 0), the expression © re- 
duces to the well-known result for a homogeneous d- 
wave BCS superconductor 0. Further, when rih = 
and V = but U ^ we find the AF state and its ex- 
pected acoustic spin cones pivoted at (tt, tt). The spin 
branches are determined by the poles of Eqn. (@J) which 
in the homogeneous case is given by the usual condi- 
tion: J7Rexo^~(q, w) — 1. This condition is no longer 
valid due to the N x N matrix structure of Eqn. Q). 
In the stripe phase without dSC, {V — 0), we show in 
Fig. |2t the imaginary part of the full spin susceptibil- 
ity. Clearly, the state contains the expected Goldstone 
modes shifted to (7r(l ± 2/iV), 7r)[l3|. The higher har- 
monics (not shown) have negligible weight compared to 
the main (7r(l ± 2/N),Tr) modes. The branches broaden 
and lose intensity as the energy is increased leaving very 
small weight near the (tt, tt) region. At low energy, the 
2D dispersion as found in constant energy scans with 
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FIG. 2: Imaginary part of the full spin susceptibility 
lmx^~{qx,qy,io) at = tt for iV = 8, [/ = 4.0t, V = 2.0t, 
and t' = 0.0 (a) and t' = -OAt (b). 



Qy ^ is dominated by circular spin cones pivoted at 
the IC points. There will be four cones when the ver- 
tical/horizontal stripe domains are equally distributed. 
There is an energy range where the cones merge and form 
weak intensity maxima at positions rotated 7r/4 from the 
Bragg IC points. In the range —0.3 < t' < 0.0, the spin 
susceptibility is similar to Fig. However, for the more 
realistic range t' < —0.3, band-structure effects cause a 
new non-Goldstone branch in the scattering response as 
shown in Fig. For t' < —OAt, this branch becomes 
fully dynamic and moves to higher energy. 

The results in Fig. 13 illustrate a general problem that 
also applies to the spin-only approaches [l4j: the outer 
branches, i.e. at \qx — 7r| > 27r/iV, are not observed 
in the low-energy NS data0, 0|. This is contrary to 
the non-superconducting La2-a;Sra;Ni04 where all four 
spin branches can be clearly seen in constant energy cuts 
through the BZ[23|- In the model calculation, the in- 
tensity along the spin branches increases monotonically 
with lowering the energy since the low-energy part of the 
damping lmxQ~ {(I, oj) qq' remain largely independent on 
the wave vector q. 

What happens when superconductivity is included? 
Then, we expect the matrix elements Xo~ i'ii^)QQ' to 
strongly depend on q at low energy. Indeed, in the ho- 
mogeneous dSC phase the particle-hole continuum con- 
tinues to w = for wave vectors connecting the nodes of 
the d-wave gap. In Fig. |3|we show representative results 
for Imx"' (q, in the coexisting phase. In the low \t'\ 
regime, the result is similar to Fig. However, for 

parameters similar to those used in Fig. ^ the result is 
very different as shown in Fig. ^jp. Clearly, at low energy 
the inner branches are considerably more intense that 
the outer. These fall within the particle-hole continuum 
and get strongly damped in agreement with the NS data 
from Christensen et This result is one of the main 

points of this paper. It is not sensitive to the size of the 
resulting gap A, but rather to the specific band-structure 
similar to the situation of a homogeneous dSC. Another 
effect induced by V is the weight around (tt, tt) which is 
significantly increased in the dSC phase. Note that the 
overall form of the spin fluctuation spectrum of Fig. Ob 
has the characteristic hourglass shape 0, 13- 

In the present approach we calculate the spin response 
in the static stripe phase relevant in the underdoped 
regime where a pseudogap is known to exist above Tc. 
Thus, the 'normal' state spin susceptibility in this dop- 
ing region is expected to be more like a thermally broad- 
ened version of Fig. as opposed to the results in Fig. 



4 




0.; 0.6 0,3 1,2 1,4 1,6 0.4 0.6 0.8 1 1,0 -.4 1,6 



(Ix [t] (Ix [tt] 

FIG. 3: lmx'^~ {qx,qy,(^) at = tt for = 8, (7 = 4.0t, 
V = 2.0t, and t' = -0.2t (a) and t' = -0.37t (b) shown 
with shghtly lower resolution than Fig. |5| Note that the 
parameters in (b) are identical to those used in Fig. 

13 In the optimally doped regime a spin gap opens due 
to the fluctuating nature of the stripes and we expect 
the intensity of the IC modes to redistribute to slightly 
above the gapQ. In the far overdoped regime the stripes 
presumably disintegrate and the picture presented here 
eventually breaks down. 

For other periodicities TV we find that at fixed U the 
spin-wave velocity is largely unchanged. Hence, the en- 
ergetic position Eres of the (tt, tt) 'resonance' is mainly 
determined by the stripe separation. This means that 
Eres increases with the doping in the underdoped regime. 
This is unlike the homogeneous dSC {U = 0, V ^ 0) 
where the intensity at (tt, tt) can be increased by band- 
structure nesting, but where Eres is solely determined by 
the maximum value of the dSC gap A which, contrary to 
experiments, decreases as the doping increases. 

Finally, note that even though the downward mode dis- 
persion at low energy (Fig. |3Jd) is qualitatively similar to 
that found in a pure dSC state, the intensity distribution 
of the IC peaks is very different [13. In the former pic- 
ture (mainly aimed at modelling YBCO and BSCCO) the 
resonance intensity is maximum at (tt, tt) and decreases 
with decreasing energy until it merges with the contin- 
uum. At lower energies the response is completely void 
due to the opening of a spin gap. 

In summary, we have calculated the electronic struc- 
ture and the associated dynamical spin susceptibility in 
the stripe phase of Eqn. For a realistic set of pa- 

rameters we find that self-consistent solutions reproduce 
salient features of widely different experimental probes 
including ARPES, STM and low-energy NS. This indi- 
cates that the stripe phase is a good starting point for 
describing the LSCO materials. In the coexisting phase 
of IC spin, charge and dSC order, the inner spin modes 
disperse toward (tt, tt) whereas the outer branches are 
strongly damped. In the future, it will be interesting to 
study in more detail the dispersion of the high-energy 
spin fluctuations within the present approach. 
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